Simple stochastic models showing strong anomalous diffusion 
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We show that strong anomalous difFusion, i.e. {\x{t)\'') ~ where qviq) is a 

nonhnear function of q, is a generic phenomenon within a class of generalized continuous- 
time random walks. For such class of systems it is possible to compute analytically z^(2n) 
where n is an integer number. The presence of strong anomalous diffusion implies that 
the data collapse of the probability density function P{x,t) = t~" F{x /t") cannot hold, 
a part (sometimes) in the limit of very small x/i", now = Yimq^ov{q). Moreover 
the comparison with previous numerical results shows that the shape of F[x/i") is not 
universal, i.e., one can have systems with the same v but different F. 
PACS number(s): 05.45. +b, 05.60.+w; 

I. INTRODUCTION 

Anomalous difFusion, i.e., when the scaling of the moments of the position x{t) is l^x^lt)) ^ t^^ with 
V > 1/2, has been observed in a rather wide class of dynamical systems, e.g., intermittent maps 
2D symplectic maps [||-^ and random velocity field&-^ as well in 2D time-dependent flow |^ and 
Hamiltonian systems (e.g. the egg-crate potential) ]lO| , |l l|] . In highly nontrivial systems, as those described 
in |^,|o|, the existence of anomalous difFusion has been established only numerically. On the other hand, 
for random shears it is possible to give an analytical criterium both for the existence of anomalous 
diffusion and for the computation of v p2[ . As far as we know, the simplest nontrivial system showing 
anomalous diffusion is the continuous-time random walk (CTRW), sometimes also called Levy walk. The 
CTRW is entirely specified by the probability density function (pdf) ^{r, r) to move a distance r in a 
time T in a single motion event. Let us assume, as in p^p^, 

^Pir,T)=P{T)P{r\T), (1) 

where P{t) is the pdf of having a flight of duration r and P{r \ r) is the conditional pdf of a displacement 
r given the flight time t. The cases corresponding to P{t) ~ r^^ and P(r \ t) = 6{\ r \ — r")/2 can be 
treated analytically |13-1^] 



If the scaling of all moments can be described by just one exponent, i.e. (a;^"(t)) ^ t^""^, a collapse of 
the pdf 's at different times is obtained exploiting the rescaling |^ 

P{x,t) ^t-'Fix/f) . (2) 

Then it also becomes clear that the value of v, in general, does not completely characterize the statistical 
properties of the diffusion process, as the function F(^) needs to be specified. 

In many cases, the use of just one exponent is not enough to describe all the moments, i.e., we have 
the so-called strong anomalous diffusion |^,|o|. For which {\x{t)\'') t'i"^'^\ where qviq) is a nonlinear 
function of q. The existence of a non-unique scaling exponent implies the failure of the data collapse for 
the pdf in the form given by Eq. 0. The best known case of a process showing strong anomalous diffusion 



is the advection of a passive scalar by a turbulent velocity field, 
that the function qv^q) is piecewise linear, i.e., 



viq q <qc 
q~c q> qc- 



In many cases M it has been observed 



(3) 



This is basically due to the existence of two mechanisms: a weak (i.e., with a unique exponent vi > 1/2) 
anomalous diffusion for the typical events, and a ballistic transport for the rare excursions (i.e., excursions 
much larger than Xtyp{t) = exp (lna;(t))). The behavior (^) suggests the vahdity of the data collapse (||) 
for the pdf core, i.e. x/f^^ not too large, and two peaks at | a; |~ t ,i.e., the footprint of ballistic events. 
The cases where strong anomalous diffusion for which ^ holds are relatively simple and surely different 
from the cases of the relative dispersion in the fully developed turbulence . 

By using only elementary techniques, in this paper we show that the bi-linear behavior for the scaling 



of the moments ( 
i.e. P(r I r) = ,5( 



3), which is present in the special case of the CTRW commonly found in the literature, 
r I — r")/2, does not hold in the general case. To show this point, we shall consider a 



generalized CTRW of the form: 



P{r 



\r) 



--S{h) 



(4) 



The inspiration for this choice comes from the multi fractal description of turbulence ||l9|] . 

The paper is organized as follows. In section II A we present the "standard" C TRW model. We then 

In section Rli 



present a simple method to find the scaling for the even order moments (section [IB), 
the model is generalized, and the same method is again applied to find the scaling. Numerical analysis 
to corroborate the analytical results of the general model are presented in Section [II together with some 
discussions related to the shape of P{x,t). Discussions and conclusions can finally be found in Section 
[V. 



II. CTRW MODELS 



Anomalous diffusion occurs when some, or all, of the hypothesis of the central limit theorem break 
down. More specifically, the system has to violate at least one of the two following conditions: 

1. Finite variance of the velocity. 

2. Fast enough decay of the auto-correlation function of the Lagrangian velocities. 

The paradigmatic model for anomalous diffusion, namely the Levy flights [ p^ violate the first condition. 
In the one dimensional case a Levy flight corresponds to the evolution in discrete time 

x{U+i) = x{U) + v,At (5) 

of the particle position x with ti^i = + At and Vi being independent stochastic variables identically 
distributed according to a Levy-stable distribution such that 

P{v) ^ for large v (6) 

where 1 < g < 3. It is easy to show that (x^) = oo for 5 < 3 and that this stochastic process shows 
anomalous diffusion being xtyp ^ ti/(9^i) > t^/^. 
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A. The "standard" CTRW model 



The existence of an infinite variance is not very pleasing from a physical point of view. This has lead 
to the introduction of the CTRW (also called Levy walks) . The idea is to relax the condition of a fixed, 
discrete time step in such a way that the process still has anomalous diffusion, but finite variance of the 
velocity. 

Firstly, we introduce the particle trajectory 

X{t) = x{t - Ti) + ViTi (7) 

where x(t) denotes the position of the particle at time t. During the random intervals t^, the particles 
move a distance r; with constant random velocity Vi independent of Ti. After each interval they choose 
new random values for r; and Vi. 

The relevant quantity to characterize the motion of the particle is the pdf ip{r, t) of having a displace- 
ment r in time r in a single motion event. This pdf is chosen in the form (|l|). In the simple case where 
Vi = ±v we have 

P(r I r) = h{\ r \ -vr) . (8) 

Taking 

P{t) - for large r (9) 

we can determine the pdf P{x,t) to be in x at time t |2^,|l[]. Actually, by introducing the probability 
density ^{x, t) to pass at location x at time t in a single motion event (and not necessarily to stop at x) 

^{x,t)=P{x\t) i dr dr ^j{r,T) ^ -6{\ X \ -vt) dr P{t) (10) 
Jt J\x\ 2 Jt 

the pdf P{x, t) can be written in the following way 

/oo j*t 
dx' dTi;{x',T)-i'{x-x',t-T) + ... (11) 
oo Jo 

the first term denotes the probability density to reach the position x at time t in a single motion event, 
the second term is the probability density to reach x at time t with one stop in x' and so on to include 
all the combinations of motion events. In the Fourier-Laplace space {x ^ k,t u) the series in Eq. ( pl| ) 
assumes the closed form 

and the behavior of {x{tY) can be calculated analytically by using the relation 



{x{tf) = 



where C ^ denotes the inverse Laplace transform [u t) [potpT] 
The results is that 
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1 < g < 2 ballistic motion 
f^-s 2 < g < 3 anomalous diffusion (13) 
t 3 < g normal diffusion. 

Thus, enhanced anomalous diffusion occurs for 2 < g < 3. For the moments of small order, which 
describe the core of the pdf, it has been shown that the asymptotic behavior gives (|a;(i)|') ^ t'^'^ with 
V = l/{g — \) [^,0. Thus the core of the pdf can be scaled as in (||) using v. The ballistic motions 
which are responsible for the different scaling of higher order moments, show up as wings on the pdf, 
which does not scale using v. 

The previous approach can be generalized [p^ to the case where 

P(r I i) = i^d r I -t"), w = ±t"-\ P(t) - with 5 > 1. (14) 

In addition one can treat more complicated situations ||l^,|l^ by considering that the particle can move 
ballistically but it can be also trapped in some structures as vortices ||2^ or chaotic islands (standard 
map or "egg-crate" potential) pl| , |2^ . 



B. Finding the scaling of the moments 

The usual method used in jl^,^,^ to find the scaling of the moments in the simple CTRW model is 
not elementary. We now present an alternative easier way to calculate the displacement moments {x'^(t)) 
and thus to characterize the anomalous diffusion. 

Let us consider a particle moving ballistically with velocity Vi during the interval times r;. The velocities 
Vi are identically distributed, independent, random variables assuming the values ±1 alternatingly. The 
intervals times are identically distributed, independent, random variables and assuming the value r 
with probability 

P(r) - T~3 with g>l and re [i,„„,r] (15) 

where tmm and T are the lowest and highest cutoffs, respectively. The reason for which we need to 
introduce such cutoffs will be clear later. The particle position at the time t can be written as 

n n 

x{t) = ^ ViTi + Vn+ie-n+i with e„+i =t-'^Ti (16) 

n being the (random) integer value for which tn < t and > t. The total time t can be rewritten as 
t = X)r=i '^i + ^n+i- Denoting by = (n) the average value of the number of time steps necessary to 
reach the time i, one has for large times t ~ N {t) and thus x{t) = X]r=i 

From simple considerations related to the symmetry of the velocity pdf under the transformation 
V I— > —V, it immediately follows that the odd-order moments {x(t)'^) are trivially zero. Conversely, 
even-order moments are nonzero and can be evaluated exploiting the following properties: (viTj) — 0, 
{viVj) (X 6ij, (TiTj) oc Sij. For times large enough, the mean squared displacement thus reads: 

(n n \ 

Similarly, for the fourth and sixth order moments the limit of large times yields 
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=N ((i;r)6) + 15Ar2 ^(„^)2^ ^(^^)4^ ^ 

15 iV^ (M')^ (18) 

and so on. 

Our attention being focused on the behavior of (x{t)'^) as a function of t, we make the substitution 
N = t/ (r) in previous expressions and (|l^) and again exploit the facts that Vi and Ti are uncorrelated 
and (^vf) = 1 to get: 



15(^)'(r^r- (19) 
For times t ^ T the leading term for (a;^"(<)) is that one proportional to i", therefore one has: 

(x(t)'^") oc I where n is integer. (20) 

V / 

which is just ordinary diffusive behaviour. Diffusive behavior in such limit is actually expected from 
general considerations. Indeed, when t ^ T the particle position at the time t can be rearranged in the 
form of a sum of almost independent displacements. If the number of the latter is large enough, central 
limit arguments apply, with the immediate consequence is that particles undergo diffusive motion. Such 
result can also be rigorously proved exploiting multiscale perturbative expansions in the (small) parameter 
T/t as done, e.g., in Ref. 

In the opposite regime, where t <^T, the system is strongly correlated, central limit arguments does not 
apply, and the final result is that non-diffusive (i.e. anomalous) regimes can occur where {x'^{t)) ^ t'^^'^i) 
with u{q) ^ 1/2. 

The possible emergence of anomalous behaviors in the limit t ^ oo can be investigated by looking at 
the dependence of moments on the cutoff T. For times shorter than T they behave as {x'^{t)) ~ t'^'^'-'^\ 

but around t ^ T the moments have a crossover to diffusive behaviour, {x'^{t)) ^ (t/ (r))''/^ (r^)''^^. 

By matching the two different regimes att = T, and using the results {t'^) ~ 7^-9+9+1 for —g + q+1 > 
and (r'^) = 0(1) for — g + g + 1 < 0, the following expressions for qv{q) are found as a function of the 
exponent g: 



(21) 



5e (1,2] 




= q 


q = 


2,4,6, 


5e (2,3] 




=q+2-g 


q = 


2,4,6, 


5e (3,4] 


qt^iq) 




q = 


2 




qt^iq) 


=q+2-g 


q = 


4,6,8, 
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and so on for higher values of g. 

From (|l]) it follows that the anomalous diffusion phase takes place for higher and higher moments as 
g increases. For q large enough one has qv{q) = (7 + constant and noting that 2i/(2) 7^ 2 one can conclude 
that ^{q) can not be constant and therefore a strong anomalous diffusion is present. 

The matching argument had been used in Ref. in the context of the multiscale method for the 
anomalous diffusion in random shear flows lal . 



C. Generalized CTRW model 

We now present a generalization of the previous model showing a strong anomalous diffusion regime 
characterized by a non-piecewise linear behavior as a function of the order q. 

As in the previous case the particle moves ballistically with random velocity Vi during the random 
interval times = i^+i — ti. We assume that has the same pdf as in (|l|) but now the velocities Vi 
assume the value ztr/* where his a random positive variable conditioned on Tj. Specifically, the conditional 
pdf P{h I r) is 

P(;i I t) oc T-^C*) (22) 

where S{h) is a positive smooth concave function. This has the effect of giving a larger variance to the 
velocity, the larger the time Tt. The function S{h) can be taylored to a special need, i.e., to mimic the 
intermittent turbulent velocity field. Here we just use a generic, simple function S{h) — /i^/(2(t^), as we 
are mainly interested in the generic properties of the model. 
The moments ((w r)^'') can now be calculated: 

{{v t)29) ^ dh dr {vT)^'^P{T)P{h I t) ~ (23) 

T p + 00 

d/ir2«''-^('') = (24) 

T 

^-9+2q+y{2q) ^ rp- g+2q+y(2q) + l ^ (^^^-^ 



where the integral / dh r^' ^'^'^^ has been evaluated exploiting the steepest descent method and we have 
defined 

y{2q) = max[2g h - S{h)] . (26) 

h 



Now using the specific shape of S{h)^ expression (g6|) takes the form: 

y{2q)^q^y[1) with y(2) = 2 a^. (27) 



Exploiting the same matching arguments as in section II B, the following expressions for the exponents 
qv{q) are obtained: 



5 6(1,2] qy(^q) = 3^a^ + q <? = 2,4,6 

ge(2,3 + y(2)] qv{q) = + q + 2 - g 9 = 2,4,6 



(28) 



5e (3 + y(2),4 + t/(4)] qv{q) ^ q/2 g = 2 

2 



qu{q)^^a-' + q + 2-g g = 4, 6, 8, 



6 



and so on, for higher values of g. 

We now see that in the general case of the CTRW, qv{q) is not just a bi-linear function, and the exact 
form depends upon the shape of S(h). The generalization to an arbitrary shape of S{h) is straightforward. 

As far as we know it is not simple to obtain the results in Eq.(pTl) and Eq.(|2^) with the method 
discussed in section II A. 



III. NUMERICAL SIMULATIONS 

We present results of numerical simulations of the general model defined by Eqs. ( p5| ) and ( |2^). The 
goal of the numerical simulations has been to verify the validity of the theoretical expectations (|21|) and 



(28) and to study the pdf of x{t) in more detail. 

In practice the pdf for the length of the jumps P(t) has to be supplemented with a lower cutoff. As it 
is the tail of the distribution which governs the scaling of the moments of x{t) the results are independent 
upon how the cutoff is made. We used a hard cutoff at r = 1, such that P(t) = for r < 1 and 

P{t) = {g- l)r"f for r > 1 . (29) 

The moments have been calculated as an ensemble average of many realizations of the process for times 
< t < T. For all the simulations presented here, T has been set to 10^. 

In Fig. |l|the scaling of the second moment (a::^(i)) with t for different values of g and a = 0.2 is shown. 
Examples have been chosen where diffusive, anomalous and ballistic behavior is expected. After an initial 
ballistic motion for short times, a transient towards scaling is taking place. For long times, clean scaling 
with an exponent corresponding to Eqs. (|2l|) is evident. 



To show how the introduction of the velocity field from section II C influences the moments, we have 
calculated the higher order moments of x{t) for a = 0,0.2 and 0.4 for g = 2.5 (Fig. ||). For the case with 
(7 = we see a clear bi-linear behavior as expected, with a cross over at Qc ~ 1.5. For low order moments 
{q < 1) the core of the pdf is important, and the behavior approches l/{g ~ 1). As is seen in the inset, 
this behavior is only to be strictly valid in the limit q 0. The scaling of the moments changes smoothly 
into the prediction 1/(5) = q + 2 — g for q > 2. Using the fact that qi'{q) is a concave increasing function 
and, for cr = 0, the slope can not be larger than 1, one obtains that the prediction given by Eq. ^ 
obtained only for even order moments, surely is valid for any moment larger than 2. For the generalized 
model (tr > 0), the fluctuations become much stronger, and it was not possible to get a clear convergence 
for large orders. However, for q = 2 the prediction in Eq. ^ and the numerical results are in perfect 
agreement. 

The characterization of the process by the scaling exponents of the moments can be seen as a way to 
probe the pdf p{x,t) of the process. Roughly speaking, the low order structure functions characterize 
the core of the pdf, while higher orders characterize the tails. In the case of "ordinary" anomalous 
diffusion, the pdfs can be rescaled according to Eq. ||, due to the fact that the whole process can be 
characterized by just one scaling exponent. However for the bi-linear or the strong diffusive cases, no 
such renormalization can be done, as more than one exponent is needed to characterize the process. As 
all the cases seem to have the same limiting behavior for g — > 0, it might be possible to make a collapse 
of the core of the pdf. In Fig. || the pdf for g = 2.5 and a = 0.2 has been rescaled using the typical value 
xtyp{t) = exp((ln(|a;(t)|)). As expected the core of the rescaled pdfs show a good overlap, while the tails 
diverge. For higher values of a the area of the core with overlap becomes smaller and smaller. 
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IV. DISCUSSIONS AND CONCLUSIONS 



By using only elementary techniques, we have shown in this paper that the strong anomalous diffusion 
appears in CTRW. Our approach, in which one computes the even moments {x{t)'^) for a system with a 
cutoff T on the pdf P(r), and then the matching of the behaviors at t > T and t <T, allows us to treat 
rather general cases. This is a relevant advantage with respect to the approach commonly used for CTRW 
which seems to us to be difficult to apply for more general cases. In the case of generalized CTRW, i.e. 
with Eq. one has a nontrivial (non bi-linear) shape of qv{q). This implies that the pdf P(a;, t) cannot 
be written in the form (||). On the other hand we found that the rescaling (H), with v — i^{0), is valid in a 
limited range of x/xtyp. It is rather natural to wonder if, at least in this limited range, the function i^(^) 
is determined by the exponent v{0). From previous works [^,|^ and the results discussed in section II], 
it seems to us that one has a negative answer: v does not determine the function F{^). 

As the shape of the function describing the scaling of the moments looks similar to what one finds 
for the relative dispersion of a passive scalar by a turbulent velocity field, one might like to use the 
the CTRW model to describe this process. There is however some differences that has to be taken into 
consideration. In fully developed turbulence, a flight of duration r can be considered as associated to 
an eddy of size I ~ r'^/^. In a process of relative diffusion, one typically has an increase of the relative 
distance with the time. This implies that a flight of duration r is preferentially followed by another flight 
with a larger duration than the previous. The easiest way for a realistic description of relative dispersion 
in fully developed turbulence should thus consist in a further generalization of the CTRW where a Markov 
ingredient, i.e. a dependence of the flight duration at the time t on the duration at the previous time, is 
introduced. It seems to us that this is a nontrivial task; a first attempt in this direction can be found in 
the recent work [G6|. 



ACKNOWLEDGMENTS 



We are particularly grateful to Paolo Muratore-Ginanneschi for very stimulating discussions and sug- 
gestions and J. Klafter for useful correspondence. The authors thank M. H. Jensen for fruitful discussions 
and all the CATS staff for the nice hospitality at Niels Bohr Institute in Copenhagen. A.M. was partially 
supported by the INFM PA GEPAIGGOl. P.C. and A.V. are partially supported by INFM (Progetto 
Ricerca Avanzata-TURBO) and by MURST (program no. 9702265437). P.C. and K.H.A. acknowledge 
the support of the European Commission's TMR Programme, Contract no.ERBFMRXCT980175, "In- 
termittency in turbulent systems" . Simulations were performed at CINECA (INFM Parallel Computing 
Initiative). 



[1] T. Geisel, J. Nierwetberg and A. Zachert, Phys. Rev. Lett., 54 (1985), 616. 

[2] J. P. Gollub and T.H. Solomon, in Proceedings of the Fritz Haber International Symposium, ed. by I. Procaccia 

(Plenum, New York, 1988) 
[3] R. Ishizaki, T. Horita, T. Kobayashi and H. Mori, Progr. of Theor. Phys., 85 (1991), 1013. 
[4] P. Leboeuf, Physica D, 116 (1998), 8. 

[5] D. Benisti and D.F. Escande, Phys. Rev. Lett. , 80 (1998), 4871. 

[6] G. Matheron and G. De Marsily, Water Resouces Res., 16 (1980), 901. 



8 



[7] M. Avellaneda and M. Vcrgassola, Phys. Rev. E, 52 (1995), 3249. 
[8] Ya.B. Zcldovich, Sov. Phys. Doki, 27 (1982), 10. 

[9] P. Castiglionc, A. Mazzino, P. Muratore-Ginanricschi and A. Vulpiani, Physica D 134, (1999), 75. 
[10] J. Klafter and G. Zumofen, Phys. Rev. E, 49 (1994), 4873. 
[11] G.M. Zaslavsky, D. Stevens, and H. Weitzener, Phys. Rev. E, 48 (1993), 1683. 
[12] L. Bifcralc, A. Crisanti, M. Vergassola and A. Vulpiani, Phys. Fluids, 7 (1995), 2725. 
[13] G. Zumofen and J. Klafter, Phys. Rev. E, 47 (1993), 851. 
[14] G. Zumofen, J. Klafter and A. Blanicn, Phys. Rev. E, 47 (1993), 2183. 
[15] G. Zumofen and J. Klafter, Europhys. Lett, 25 (1994), 565. 
[16] M. Schlesinger, B.J. West and J. Klafter, Phys. Rev. Lett, 58 (1987), 1100. 
[17] G. Boffotta, A. Cclani, A. Crisanti and A. Vulpiani, Europhys. Lett. 42, (1999), 177. 

[18] E. MontroU and M. Schlesinger, in Studies in Statistical Mechanics edited by E.W MontroU and J.L. Lcbowitz, 

vol. 11, p. 1, North- Holland, Amsterdam, 1984. 
[19] T. Bohr, M.H. Jensen, G. Paladin and A. Vulpiani Dynamical Systems Approach to Turbulence (Cambridge 

University Press, Cambridge UK 1998) 
[20] G. Zumofen, J. Klafter and M.F. Schlesinger, in Anomalous diffusion. From Basics to Applications edited by 

R.Kutner, A.Pekalski and K. Szanajd-Weron, Lecture Notes in Physics, Vol. 519, p. 15, Springer (1999). 
[21] G. Zumofen, J. Klafter and M.F. Schlesinger, in Ghaos-The interplay between stochastic and deterministic 

behaviour editeA by P. Garbaczewski, M. Wolf and A. Weron, p. 183, Springer, (1996). 
[22] T.H. Solomon, E.R. Weeks and H.L. Swinney, Phys. Rev. Lett, 71 (1993), 3975. 
[23] T. Geiscl, A. Zacherl and G. Radons, Z. Phys., B71 (1988), 117. 
[24] J.P. Bouchaud and A. Georges, Phys. Rep., 195 (1990), 127. 

[25] J.-P. Bouchaud, A. Georges, J. Koplik, A. Provat and S. Redner, Phys. Rev. Lett, 64 (1990), 2503. 
[26] I.M. Sokolov, A. Blumen and J. Klafter Europhys. Lett, 47 (1999), 152 

1e+10 
1e+09 
1e+08 
1e+07 
1e+06 
A 100000 
V 10000 
1000 
100 
10 

1 

0.1 

0.1 1 10 100 1000 10000 

t 




9 



FIG. 1. The scaling of the displacement (^x^{t)^ for three different values of g, each corresponding to diffusive 
behavior {g = 4.0), anomalous diffusion {g = 2.5) and ballistic motion {g = 1.5). The value of a is 0.2 for all the 

cases. In the inset the same is shown, but this time roscalcd with the expected behavior (x^) oc t^'^^'^\ such that 
a scaling in accordance with the theoretical prediction corresponds to a horizontal line. 
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FIG. 3. The pdf of x{€) at three difTerent times rescaled according to Eq. (Q) for a situation with g — 2.5 and 
CT = 0.2 corresponding to the anomalous regime. The pdf have been scaled with the typical value, Xtyp{t). 



11 



